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The charm content in r\ is estimated using the QCD U(l) axial anomaly applied to 

" ~~ * s 99 —* w s. We find / , f» —3.1 MeV for m c = 1.3 GeV. Our estimate agrees 
n 

with other independent methods by Feldmann and Kroll, and Araki, Musakhanov 
and Toki. The resulting branching ratios for B — > n K in the generalized factor- 
ization model are marginally consistent with the CLEO data, although they are 
lower than the data by ~ 2. 

1 Introduction 

In the year of 1997, the CLEO collaboration reported measurements in a num- 
ber of exclusive two-body non-leptonic decays of the type B =-> iji/12, where h\ 
and hi are light mesons and the inclusive decay B^ — > rj'X s . an In particular, 
large branching ratios into the final states including rf are reported : u 

BR{B ± -► ?y' + X s ) = (6.2 ± 1.6 ± 1.3) x 1(T 4 (1) 

(for 2.0 GeV < p r < < 2.7 GeV), 

BR(B ± rf + X ± ) = (7.1^1 ± 1.0) x 1CT 5 , (2) 

BR(B° -» rj'K ) = (5.3±1;1 ± 1.2) x 10" 5 . (3) 

They did not observe any decay involving the (i]K) or (77, 7;')^*_modes. These 
measurements have stimulated a lot of theoretical activity, Utli both in the 
inclusive and exclusive decays involving 77 mesons. 

Such unexpectedly large branching ratios of B decays into the final states 
with an 77 meson led to an idea that the charm content in 77 meson might be 
very large. The amplitude for b — > s(cc) — > s(j]',i]) can be parametrized as 
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(??'|c7p75c|0) = —ifZq.ni and similarly for the 77 case. The relevant quantity 
/ c (/) is often referred to as the charm content of the rf-'K This quantity is a 
priori not known, but it was claimed that the data (l)-(3) can be accomodated 
if fZ is as large as (50 — 180) MeV. Qo However, this is uncomfortably large 
compared to the typical meson decay constants, e.g., f n = 132 MeV. Later 
it turned out that they could be determined in a number of different, ways, 
also including the B-decays combined with the QCD E/a(1) anomalylla being 
discussed here. 

In this talk, I'd like to present another method for computing the con- 
tribution of the amplitudes b — > s(gg) — > s(r}',rj). This method is based on 
calculating the amplitude for the chromomagnetic penguin process b — > sgg, 
followed by the transitions gg — > (77', 77) which are calculated using the QCD 
anomaly, determining both the sign and magnitude of these contributions. The 
branching ratios for B ± -> (77', 77) (A"*, K*±) and B° -> (77', rj)(K°, K*°) based 
on the QCD-anomaly method are calculated in this letter and compared with 
the present CLEO measurements and with the ones in Ref. [5]. We find that 
the theoretical branching ratios for B ± — ► rj'K^ and B° — > r/'K are almost 
equal and both are in the range (2 — 4) x 10 -5 , in agreement with the estimates 
in Ref. [5]. 

2 Estimate of b — > (r],r]')s via QCD anomaly 

Let us consider the charm- anticharm pair into two gluons, followed by the 
transition gg — » 7/ ' (see Fig. 1). The first part of this two-step process, i.e. 
b — > s(cc — > g{k\)g{k2)) which amounts to calculating the charm-quark-loop 
m which two gluons are emitted, has been worked out by Simma and Wyler 
in the context of a calculation in the full theory. Their result is readily 
translated to our effective theory approach and can be compactly written as a 
new (induced) effective Hamiltonian H^, 

xGf(DpG afi ) a s 7 "(l- 75 )6 , (4) 

with G^v — ^e^afjG 01 ^ (£0123 = +1)< The are modified Wilson coeffi- 
cients for the current-current operators as described in Ref. [5] . In this formula, 
which holds for on-shell gluons (q 2 — (k\ + k^) 2 = 2ki ■ k2), the sum over color 
indices is understood. The function Ais(q 2 /m 2 .) is defined as 

2 

, for < z < 4 . (5) 
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z 



n — 2 arctan( 1) 

2 



1/2 



2 



Note that the b — ► sgg calculation brings in an explicit factor of a s . However, 
this explicit a s factor gets absorbed into the matrix element of the operator 
resulting from the anomaly. So, to the order that we are working, we use the 
coefficients Cf and Cf in Eq. (@). The uu contribution in Fig. 1 is suppressed 
due to the unfavourable CKM factors. The ti contribution is included in the 
effective Hamiltonian via the bsgg piece present in the operator O g . However, 
in the factorization framework, the bsgg term in O g does not contribute to 
the decays discussed. So, the cc contribution in Fig. 1 is the only one that 
survives. 

Working out the hadronic matrix element of Eq. (|]) using factorization, 
we now need to evaluate the matrix elements : 

(^\Gf(D p G a ^ a )\0), (6) 

which can be written as 

G" 13 (DpGa^a — dp{G^G ail<a ) — {DpG al3 ) a G ail ^ a . (7) 

Now we can discard the second term since it is suppressed by an additional 
power of g s which follows on using the equation of motion, and furthermore, 
the first term is enhanced by N c in the large N c limit. The matrix elements of 
) are related to those of GG : 

^(/)|GfG^ )O |0) = ^(/)|G^G^, a |0) . (8) 

The conversion of the gluons into r\ and rf is described by an amplitude which 
is fixed by the SU(3) symmetry and the axial LL/1) current triangle anomaly. 
The matrix elements for GG can be written asEJ 

( v O\^GfG a p, a \0) = rnl ( , ) f% e) . (9) 
In Eq. (||) the decay constants and read 

/«= ^| cos 8 - ^| sin O , # = ^= sin ^ + ^ cos O , (10) 

where the coupling constants /g, fo and the mixing angles 08 and 0o have 
been introduced earlier. We follow here the two-angle (77,77') mixing formal- 
ism of Ref. [21], where the mass eigenstates \rj) and I77') have the following 
decompositions: 

I?;) = cos0 8 |?78) - sin0 o |77o), 

1 77') = sin0 8 |f7 8 ) + cos0 o |?7o)- (11) 
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Collecting the individual steps, the matrix elements in Eq. (6) can be written 
as 

9 

771 

(v ( 'Hq)\^Gf(D p G a , a )\0)= iq ,^ff; (l> . (12) 

One would have naively expected that the gluonic matrix elements are small 
since they contain an extra factor of a s . However, as shown by Eq. (p"2|), 
this is obviously not the case and the gluon operator with a s as a whole is 
responsible for the invariant mass of the rf 1 ' mesons. Also, the combination 
entering in Eq. ( |l2] ) involving the product of a s and the gluon field operators 
is independent of the renormalization scale. 

Before closing this section, let me comment on the alternative method to 
estimate f c , using the i] c — 77 mixing (Fig. [2] (a)) : 

(13) 

For example, in Ref. [5], these quantities were determined from the decays 
J/ip — > (77, tj' ,T)c)~fi extending the usual (77, ?/)-mixing formalism E3 to the 
(77c 77', 77) system. Using the measured decay widths for the decays J/tp -* 
(v,v',Vch and (r^, 7/, t?) - 77 yields ~ 5.8 MeV and \f<\ ~ 2.3 MeV. U 
However, in such a scheme, the intermediate 77 c is far off-shell, and one has to 
know the form factor for J/tp — r\ c ~ 7 vertex. Moreover, this picture is not 
consistent with electromagnetic gauge invariance, since at the quark level, the 
above transition occurs through cc — * 7 + gg followed by gg — ► 77 (Fig. || (b)). 
Note that there is a diagram in which the photon is emitted in the middle of 
two gluons. Such diagram cannot be represented by the simple rj c — rj mixing 
picture. One cannot simply neglect this diagram either, since it would violate 
electromagnetic gauge invariance. Therefore, the usual picture of ?7 C — 77 mix- 
ing for J/tJj — > 77 7 may not be a good one, although it is widely used in the 
literature. 



3 Estimate of the decay rate for B — > (77', rj)(K, K*) 

In order to compute the complete amplitude for the exclusive decays, one has to 
combine the contribution from the decay b — > s(cc) — * s(gg) — > sr]( ' discussed 
in the previous section with all the others arising from the four-quark and 
chromomagnetic operators, as detailed in Ref. [5]. The resulting amplitudes 
in the factorization approximation are listed in Ref. [18]. By comparing two 
expressions in Refs. [5] and [18], we can make the following identification : 
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Therefore, we have a simple relation between the decay constants 
introduced in the intrinsic charm content method, and the form factor Ai$ 
entering via the operator in Eq. (Q) . The idea of intrinsic charm quark content 
of 7/ and 77 and the contribution of the operator in Eq. (|J) are related since 
this operator comes from the charm quark loop. Using the best-fit values 
of the (77, ?7')-mixing parameters from Ref. [23], yielding 6*g = —22.2°, 80 = 
-9.1°, / 8 = 168 MeV,/ = 157 MeV, which in turn yields = 63.6 MeV 

and = 77.8 MeV, the relations in Eq. (|4j) give /° 3.1 MeV (-2.3 

MeV) and 1.2 MeV (-0.9 MeV), with m c having the value 1.3 GeV 

(1.5 GeV). In our approach the relative signs of the contributions from 0\ and 
O2 to the other contributions are determined; we obtain the negative- f°, (and 
/°) solution of the two possible ones. Our estimate-pf f v > is consistent with 
the results obtained by Feldmann, Kroll and Stech,c3 

f v > = -(6.3 ± 0.6) MeV, f n = -(2.4 ± 0.2) MeV, (15) 

which is based on a new scheme for the flavor mixing and the axial anomaly. 
Also there appeared Ref. [25] by Araki et at, who did an independent calcu- 
lation in the same way as Haperin and Zhitnitsky and got 

/ ( , c) = -(12.3 - 18.4) MeV, (16) 

after correcting some mistakes in the orginal calculations in Ref. [7] . By now 

r (c) 

it is amusing that three independent methods for estimating / , agree on its 
sign and the size within a factor of 2 — 3. So one can say that the issue of the 
charm content in 77 is now settled with reasonable confidence. 

We plot the resulting branching ratios BR(B^ — > r/K a± ) and BR(B° — > 
rj'K ) as functions of the parameter £ in Figs. 3 and 4 for three-different sets 
of CKM elements p and 77 in the Wolfenstein parametrization : 

0,?y) = (0.05,0.36), (0.30,0.42), and (0,0.22) (17) 

The branching ratios for the neutral i?-meson decays are averages with the 
corresponding charged conjugated decays in the figures. The numerical results 
for the branching ratios depend on various input parameters, the details of 
which can be found in Ref. [18]. Let us discuss the sensitivity on various 
parameters in brief. The branching ratios show a mild dependence (of order 
10%) on the CKM parameters. However the branching ratios depend on the s- 
quark mass very sensitively, because the amplitude for B — > 77 K contains terms 
with the l/m s factor. In this work, we used m s (2.5 GeV) = 122 MeV. However, 
if the present high values of the branching ratios for B ^ — > rj'K ^ continue 
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to persist, one might have to consider smaller values of mj. We also expect 
progress in calculating quark masses on the lattice, sharpening the theoretical 
estimates presented here. Dependence on £ amounts to between 20% and 35% 
depending on the other parameters if one varies £ in the range < £ < 0.5. 
In all cases, the branching ratios are larger for £ = 0. Taking into account the 
parametric dependences just discussed, we note that the theoretical branching 
ratio BR(B ^ — > rj'K ^) are uncertain by a factor 2. For the ratio of 
the branching ratios BR(B ± -> rj'K ± )/BR(B° -> rfK ), which is useful 
quantity since it is practically independent of the form factors and most input 
parameters, the residual uncertainty is due to the CKM-parameter dependence 
of this ratio. It is estimated as about 10%. We get (for < £ < 0.5) 

BR(B ± -> r/K*) , N 

- =0 - 9 - L02 - ( 18 ) 

The present experimental value of this ratio as calculated by adding the exper- 
imental errors in the numerator and denominator in quadrature is 1.34 ±0.85. 
One can also study other modes involving K* instead of K. Summarizing the 
results of Ref. [18], the branching ratios for the decay modes B^ — > r/K^ and 
B° — > r/K are smaller compared to their ^'-counterparts by at least an order 
of magnitude, namely, BR(B^ — ► r/K ± ) = (1 — 2) x 10~ 6 and a similar value 
for the neutral B decay mode. On the other hand, the branching ratios for the 
decay modes B* 1 — > ■q(K ± ,K* ± ) and B° — > r](K°,K*°) are all comparable to 
each other somewhere in the range (1 — 3) x 10 -6 . 



4 Concluding Remarks 

In this talk, we have presented an independent estimate of the charm con- 
tent of rj( ' meson using the process b — * s(cc) — * s(gg) — * 5(77', 77) and QCD 
anomaly. We could fix its sign as well as the size. Our results are in good 
agreement with other results based on independent methods, and thus we can 
be confident by now that the charm content in rf •* meson is indeed small. One 
could further study the branching ratios in B — » (77', r])(K, K*) in the general- 
ized factorization approximation. Our result is marginally consistent with the 
CLEO data, considering various theoretical uncertainties in this kind of game 
(Figs. 3 and 4). Also our result is drastically different from the ones which 
follow in other scenarios. Hence, ongoing and future experiments will be able 
to test the predictions of the present approach as well as of the competing ones, 
such as models based on the dominance of the intrinsic charm contributions in 
77', as suggested in Refs. [7,8], or models in which damipant role is attributed 
to the soft-gluon- fusion process to form an 77 or 77'. 113-0 
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Figure 1: Fcynman diagram contributing to the processes b — > s(cc) — > s(gg) — > s?/ ' in 
the full and effective theory. The lower vertex in the diagram on the right is calculated with 
the insertion of the operators 0\ and 0| in the effective Hamiltonian approach; the upper 
vertex in both the full and effective theory is determined by the QCD triangle anomaly. 
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Figure 2: Fcynman diagram contributing to the processes J/tp —* 77J (a) in the »7 C — eta 
mixing picture, and (b) in the parton picture. 



7 




Figure 3: The branching ratio BR(B^ — > rj K^) plotted against the parameter £. The 
lower three curves correspond to the value m s (2.5 GeV) = 122 MeV and the three choices 
of the CKM parameters: p = 0.05, T) = 0.36 (solid curve); p = 0.30, r) = 0.42 (dashed 
curve); p = 0, r\ = 0.22 (dashed-dotted curve). The upper two curves correspond to the 
value mj(2.5 GeV) = 100 MeV, p = 0.05, r) = 0.36 and £= , = -2.3 MeV from the QCD- 
anomaly method (dotted curve) and f c , = —5.8 MeV from|3 (long-short dashed curve). The 
horizontal thick solid lines represent the present CLEO measurements (with ±ltr errors). 
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Figure 4: The branching ratio BB.(B° — > r/ K°) plotted against the parameter £• The 
legends are the same as in Fig. 3, and the horizontal thick solid lines represent the present 
CLEO measurements (with ±ltr errors). 
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